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ABSTRACT 

UGC 4483 is a nearby Blue Compact Dwarf (BCD) galaxy HST observations have resolved the galaxy into single stars and this has led 
to the derivation of its star formation history and to a direct estimate of its stellar mass. We have analysed archival VLA observations 
of the 21 -cm line and found that UGC 4483 has a steeply-rising rotation curve which flattens in the outer parts at a velocity of 
~20 km s _1 . Radial motions of ~5 km s _1 may also be present. As far as we know, UGC 4483 is the lowest-mass galaxy with a 
differentially rotating H I disk. The steep rise of the rotation curve indicates that there is a strong central concentration of mass. We 
have built mass models using the HST information on the stellar mass to break the disk-halo degeneracy: old stars contribute ~50% of 
the observed rotation velocity at 2.2 disk scale-lengths. Baryons (gas and stars) constitute an important fraction of the total dynamical 
mass. These are striking differences with respect to typical dwarf irregular galaxies (dlrrs), which usually have slowly-rising rotation 
curves and are thought to be entirely dominated by dark matter. BCDs appear to be different from non- starbursting dlrrs in terms of 
their H I and stellar distributions and their internal dynamics. To their high central surface brightnesses and high central H I densities 
correspond strong central rotation- velocity gradients. This implies that the starburst is closely related with the gravitational potential 
and the concentration of gas. We discuss the implications of our results on the properties of the progenitors/descendants of BCDs. 

Key words, dark matter - galaxies: individual: UGC 4483 - galaxies: dwarf - galaxies: starburst - galaxies: kinematics and dynamics 
- galaxies: evolution 



1. Introduction 

The mechanisms that trigger strong bursts of star formation in 
galaxies are poorly understood. In the Local Universe, starburst 
activity is mostly observed in low-mass galaxies, which are usu- 
ally classified as blue compact dwarfs (BCDs) (e.g. Gil de Paz 
et al. 2003), amorphous dwarfs (e.g. Gallagher & Hunter 1987), 
or HII galaxies (e.g. Taylor et al. 1995). Hereafter, we will refer 
to any starbursting dwarf galaxy as a BCD. Several studies (e.g. 
Gil de Paz & Madore 2005, Tosi 2009 and references therein) 
have shown that BCDs are not young galaxies undergoing their 
first burst of star formation (as suggested by Searle & Sargent 
1972), as they also contain old stellar populations with ages >2- 
3 Gyr. In particular, HST has made it possible to resolve nearby 
BCDs into single stars and to derive colour-magnitude diagrams 
deep enough to provide the following information: i) accurate 
distances of the galaxies, ii) a direct estimate of their total stel- 
lar mass, and iii) their star formation history (SFH) (e.g. Tosi 
2009). These SFHs show that the starburst is a short-lived phe- 
nomenon, typically sustained for a few 10 8 yr (McQuinn et al. 
2010a). Thus, BCDs are transition-type dwarfs but the nature of 
their progenitors and descendants remains unclear. In particular, 
it is not known whether there are evolutionary connections with 
dwarf irregulars (dlrrs), spheroidals (dSphs), and/or ellipticals 
(dEs) (e.g. Papaderos et al. 1996; van Zee et al. 2001). 

There are striking differences between BCDs and other types 
of dwarf galaxies: i) the old stellar component of BCDs gene- 
rally has a smaller scale-length and higher central surface bright- 
ness than dlrrs and dEs/dSphs (e.g. Papaderos et al. 1996; Gil de 
Paz & Madore 2005); ii) BCDs have strong concentrations of 
H I within the starburst region, where the column densities are 
typically 2-3 times higher than in dlrrs (e.g. van Zee et al. 1998, 



2001); iii) BCDs have steep central velocity gradients that are 
not observed in dlrrs (e.g. van Zee et al. 1998, 2001). The steep 
velocity gradients may signify a steeply-rising rotation curve 
(van Zee et al. 2001; Lelli et al. 2012), high velocity disper- 
sion, or non-circular motions (e.g. Elson et al. 2011b). Detailed 
studies of the gas kinematics are needed to determine the inner 
shape of the rotation curve. Recently, Lelli et al. (2012) studied 
the BCD prototype I Zw 18 and found that it has a flat rotation 
curve with a steep rise in the inner parts, indicating that there is 
a high central concentration of mass. Such a mass concentration 
is not observed in typical dlrrs. This points to a close connec- 
tion between the starburst and the gravitational potential. It is 
also clear that a BCD like I Zw 18 cannot evolve into a typical 
dlrrs at the end of the starburst, unless the central concentra- 
tion of mass is removed. It is important, therefore, to investigate 
whether all BCDs have steeply-rising rotation curves, and to de- 
termine the relative contributions of gas, stars and dark matter to 
the gravitational potential. 

We present a detailed study of the gas kinematics of 
UGC 4483, a starbursting dwarf galaxy located in the M81 
group and resolved into individual stars by HST (Dolphin et al. 
2001; Izotov & Thuan 2002). UGC 4483 is extremely metal- 
poor (12+log(0/H)^7.5, see Skillman et al. 1994 and van Zee & 
Haynes 2006) and may be classified as a "cometary" BCD, as its 
high-surface-brightness starburst region is located at the edge of 
an elongated low-surface-brightness stellar body (see Fig. 1, top- 
left). Previous HI studies (Lo et al. 1993; van Zee et al. 1998) 
showed that the galaxy has an extended H I disk, with a strong 
H I concentration near the starburst region and a steep central 
velocity gradient, but the H I kinematics was not studied in de- 
tail. We analysed archival VLA data and were able to derive a 
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Table 1: Properties of UGC 4483. 



a (J2000) 


08 h 37 m 3. s l +0. s 5 


6 (J2000) 


69°46'31" + 2" 


Distance (Mpc) 


3.2 + 0.2 


V sys (km s" 1 ) 


158 + 2 


Position Angle (°) 


U + J 


Inclination Angle (°) 


58 + 3 


V rot (km s" 1 ) 


19 + 2 


M dyn (10 7 M ) 


16 + 3 


M, (10 7 M ) 


1.0 + 0.3 


M H i (10 7 M Q ) 


2.5 + 0.3 


L B (10 7 L Q ) 


1.4 


L R (10 7 L ) 


0.9 



Notes. Luminosities were calculated using the apparent magnitudes 
from Gil de Paz et al. (2003), the distance from Dolphin et al. (2001) and 
the solar absolute magnitudes from Binney & Merrifield (1998). The 
stellar mass was calculated integrating the SFH from McQuinn et al. 
(2010a) and assuming a gas recycling efficiency of 30%. The dynami- 
cal mass was calculated taking into account the pressure-support. 



rotation curve which we used to investigate the distributions of 
luminous and dark matter in this galaxy. 

2. Data reduction & analysis 

We analysed H I data taken from the VLA archive. The obser- 
vations have been carried out with the B and C arrays, and are 
described in van Zee et al. (1998). The correlator was used in 
2 AD mode, with a total bandwidth of 1.56 MHz (-330 km s" 1 ). 
An on-line Hanning taper was applied to the data, producing 127 
spectral line channels with a width of 12.3 kHz (~2.6 km s" 1 ). 

We interactively flagged, calibrated, and combined the raw 
UV data using the AIPS package and following standard VLA 
procedures. The UV data were mapped using a robust weighting 
technique (Briggs 1995) and a Gaussian baseline taper to atten- 
uate the longest baselines. After various trials, we chose a robust 
parameter of -1 and a taper FWHM of 40 kA; these parameters 
minimize sidelobes and wings in the beam profile and lead to a 
datacube with an angular resolution of 5". 7 x 4". 5. 

After the Fourier transform, the data analysis was contin- 
ued using the Groningen Imaging Processing SYstem (GIPSY) 
(van der Hulst et al. 1992). A continuum map was constructed 
by averaging line-free channels, and subtracted from the data- 
cube. The channel maps were cleaned (Hogbom 1974) down 
to 0.3 cr, using a mask to define the search areas for the clean- 
components that were then restored with a Gaussian beam of the 
same FWHM as the antenna pattern. The mask was constructed 
by smoothing the datacube both in velocity (to 10.4 km s" 1 ) and 
spatially (to 20") and clipping at ~3cr s (where cr s is the rms 
noise in the smoothed cube). In order to improve the S/N, the 
cleaned datacube was smoothed in velocity to a resolution of 
5.2 km s" 1 and spatially to 10", providing a 3<x column density 
sensitivity of 6xl0 19 atoms cm" 2 per 2.6 km s _1 -wide channel. 

A total H I map was constructed by summing the signal in- 
side the clean-mask; a pseudo-3<x contour was calculated fol- 
lowing Verheijen & Sancisi (2001). A velocity field was derived 
by fitting a Gaussian function to the HI line profiles. Fitted 
Gaussians with a peak intensity less than 3cr and a FWHM 
smaller than 5.2 km s" 1 were discarded. The HI line profiles 
are quite broad and asymmetric, thus the velocity field derived 
from the Gaussian fitting provides only an overall description 
of the galaxy kinematics. Our kinematical analysis is based on 



three-dimensional (3D) models of the observations (Sect. 3.2) 
and not merely on the two-dimensional (2D) velocity field. 

3. Results 

3.1. Hi distribution and kinematics 

Figure 1 shows the total H I map of UGC 4483 at a resolution of 
10" (top-right) and a R-band image at the same scale (top-left). 
The H I distribution is lopsided and closely resembles the optical 
morphology. There is a strong H I concentration near the star- 
burst region to the North. At 5. "7x4". 5 resolution (-80 pc), the 
peak column densities are ~5xl0 21 atoms cm -2 (-40 M© pc -2 ) 
(same as found by van Zee et al. 1998). 

The H I kinematics of UGC 4483 is illustrated in Fig. 1 (bot- 
tom) and Fig. 2 (top). The HI emission shows a velocity gradi- 
ent along a position angle P.A.~0°, which roughly corresponds 
to the optical major axis of the galaxy, suggesting that there is 
a rotating H I disk. The velocity field indicates large-scale diffe- 
rential rotation. It also shows, however, large-scale asymmetries. 
In the next section, we derive the rotation curve of UGC 4483 
and discuss possible non-circular motions. 

3.2. Rotation curve and non-circular motions 

Rotation curves of disk galaxies are usually derived by fitting a 
tilted-ring model to a velocity field (e.g. Begeman 1987). In the 
case of UGC 4483, there are severe limitations because of the 
asymmetries in the velocity field. Van Zee et al. (1998) modeled 
the velocity field of UGC 4483 and obtained a rough estimate 
of the dynamical mass. We derived a rotation curve using the 
following approach. As a first step, we fit a tilted-ring model 
to the velocity field to obtain an initial estimate of the rotation 
curve. Then, this rotation curve was used as input to build a 3D 
kinematic model and subsequently corrected by trial and error to 
produce a model-cube that matches the observations. 

When fitting a tilted-ring model to the velocity field, we used 
a ring width of 10" (1 beam), thus the points of the rotation 
curve are nearly independent. The points of the velocity field 
were weighted by cos 2 (6), where is the azimuthal angle in the 
plane of the galaxy. We kept the centre (xo,yo) and the incli- 
nation i fixed, using the values derived from the optical image 
by fitting ellipses to the outermost isophotes (see table 1). An 
inclination of 58° is formally a lower limit, because the stellar 
component may be thick; if the disk is assumed to be 5° more 
edge-on, the rotation velocities would decrease by only -5%. 
We also fixed the PA. assuming the value of 0° suggested by the 
HI morphology (Fig. 1, top), which is consistent with the op- 
tical PA. (^ -10°) within the uncertainties. As a first step, we 
determined the systemic velocity V sys , by taking the mean value 
over all the rings. Then, we fixed V sys and determined the rota- 
tion velocity V rot and the radial velocity V ra d for each ring. The 
rotation curve has an amplitude of -18-20 km s" 1 (Fig. 4), while 
the radial velocities are ~4-5 km s" 1 . 

As a final step in the derivation of the rotation curve, we built 
3D kinematic models, similarly to Swaters et al. (2009) and Lelli 
et al. (2010, 2012). The disk kinematics is assumed to be axi- 
symmetric, while the H I distribution is clumpy, i.e. the surface 
density varies with position as in the observed H I map. The pro- 
cedure is as follows: i) a disk with uniform surface density and 
constant thickness is constructed by fixing the velocity disper- 
sion <thi , the rotation velocity V rot and the radial velocity V ra d at 
every radius; ii) the disk is projected on the sky using the geo- 
metrical parameters (xq, yo), V sys , PA., and i and a model-cube 
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Fig. 1: Top-Left: R-band image (from Gil de Paz et al. 2003). Top-Right: total HI map at 10" resolution. The contour values are 1.7, 
3.4, 6.8, 13.6, 27.2 x 10 20 atoms cm" 2 . Bottom: channel maps at 10" resolution. Contours are at 1.5, 3, 6, 12, 24 x <x, where cr=0.66 
mJy/beam (~2xl0 19 atoms cm -2 ). In each panel, the cross and the circle show the galaxy centre and the beam size, respectively. 



is created; iii) the model-cube is convolved with the observa- 
tional beam; and iv) the H I line profiles are rescaled to repro- 
duce the observed H I map, i.e. the flux density is recovered at 
every spatial pixel. For the geometrical parameters, we used the 
values in table 1. For the vertical distribution, we assumed an 
exponential-law exp(-z/zo) with zo=100 pc. We also assumed 
that <thi =8 km s" 1 over the entire disk; a mean velocity disper- 
sion higher than 10 km s" 1 can be ruled out by comparing the 
shape of position- velocity (PV) diagrams obtained from the mo- 
dels and the observations. The values of zo and <thi are slightly 
degenerate, but do not significantly affect the final results. For 
the rotation curve, we used the rotation velocities obtained by fit- 
ting the velocity field, then we corrected them by trial and error 
to obtain a model-cube that matches the observations. We find a 
good match by increasing the first four points by ~3 km s" 1 . This 
rotation curve is shown in Fig. 4. The comparison between mo- 



dels and observations has been done using velocity fields (see 
Fig. 2) and PV-diagrams (see Fig. 3), Note that the model ve- 
locity fields were derived from the model-cubes and, therefore, 
include the effects of spatial and spectral resolution, velocity dis- 
persion, and possible non-circular motions. Because in our mo- 
dels we initially assumed a uniform H I disk, some aspects of 
the effects of beam- smearing are not fully included if the H I 
distribution varies rapidly on small scales. However, we used 
PV-diagrams at the full resolution (5. "7x4". 5) and compared 
them with those at 10" to check that beam- smearing effects do 
not significantly affect our rotation curve. 

Figure 4 (filled-circles) shows that UGC 4483 has a steeply- 
rising rotation curve that flattens in the outer parts at a velocity 
of -18-20 km s _1 . This is a striking difference with respect to 
dlrrs of similar mass, as they usually have a slowly-rising rota- 
tion curve close to a solid body. As far as we know, UGC 4483 is 



3 



Lelli et al.: Dynamics of starbursting dwarf galaxies II: UGC 4483. 



Observations 




R.A. (J2000) R.A. (J2000) 

Fig. 2: Velocity fields derived from different 3D kinematic models and the observations. The rotation curve used to build the models 
is shown in the inset on the left (x-axis: radius in arcsec, y-axis: rotation velocity in km s" 1 ). Contours range from 140 to 181.6 
km s" 1 , with steps of 5.2 km s" 1 . The circles show the beam size (10"). See Sect. 3.1 for details. 
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SOLID BODY DIFF.ROT. DIFF.ROT. + RAD. OBSERVATIONS 
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Fig. 3: Position- velocity diagrams derived for different 3D kinematic models and from the observations. The slices are taken along 
the major and minor axes. Contours are at 1.5, 3, 6, 12 x <x, where <x=0.66 mJy/beam. Squares show the rotation curve used to build 
the models. See Sect. 3.1 for details. 
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Fig. 4: Observed HI rotation curve (blue-filled circles) and 
asymmetric-drift-corrected rotation curve (red-open circles). 



the lowest-mass galaxy with a differentially rotating H I disk. To 
further test the validity of this result, we built 3D kinematic mo- 
dels assuming a solid-body rotation curve and compared them 
with the differentially-rotating disk-model and with the obser- 
vations. The results are shown in Fig. 2 and Fig. 3. It is clear 
that a flat rotation curve provides a better match of the obser- 
vations than a solid-body one, as the iso-velocity contours of 
the observed velocity field display a curvature that is typical 



of differential rotation. The same conclusions are reached by 
comparing PV-diagrams obtained from the models and the ob- 
servations along the major axis (Fig. 3, top). The observed PV- 
diagram shows a "flattening" on the southern-approaching side 
of the disk, that is reproduced only by a differentially rotating 
disk. On the northern-receding side, the "flattening" is less clear 
as the disk is less extended, but a differentially rotating disk is 
still preferable than a solid body one. 

A simple rotating disk, however, cannot reproduce the asym- 
metries present in the observed velocity field. In particular, the 
kinematic minor axis, defined by the contours close to the sys- 
temic velocity, is not orthogonal to the major axis. This may be 
due either to radial motions (e.g. Fraternali et al. 2002) or to an 
oval distortion of the gravitational potential (e.g. Bosma 1978). 
These asymmetries can be seen also in the PV-diagram taken 
along the geometrical minor axis, i.e. in the direction perpendi- 
cular to the major axis (Fig 3, bottom). Thus, we improved our 
models by adding a constant radial component of ~5 km s" 1 , 
which is also indicated by the tilted-ring fit to the velocity field. 
A model with solid-body rotation plus radial motions is still not 
acceptable (Fig. 2, bottom-right), whereas a model with differen- 
tial rotation plus radial motions reproduces most of the features 
in the data (Fig. 2, middle-right and Fig. 3, bottom). 

The errors on the rotation curve have been estimated as 
^"rot ~ °"fit + ^"asym' where cr^ is the formal error given by the 
tilted-ring fit and <r^ sym is an additional uncertainty due to the 
asymmetries between the approaching and receding sides, that is 
estimated as <T aS ym = (V ro t,ap P - V rot ,rec)/4 (Swaters et al. 2009). 
We point out that the observed H I line profiles are quite broad 
and, therefore, it is very difficult to trace the rotation curve in 
the innermost parts. We compared models that have an inner 
solid-body rise and flatten at different radii and estimated that 
the rotation curve must flatten between -0.3 and -0.6 kpc, giv- 
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Fig. 5: Left: Velocity fields derived from disk-models with diffe- 
rent geometrical parameters: optical centre (cross), P.A.=0°, and 
radial motions of ~5 km s" 1 (top panel, same as Fig. 2 middle- 
right); shifted dynamical centre (circle), P.A.=10°, and no ra- 
dial motions (bottom panel). Contours are the same as in Fig. 2. 
Right: differences between the observed and the model velocity 
fields. Contours are at ±2, ±4, ±6, ±8 km s" 1 . The filled circle 
shows the beam. 



ing an inner rotation- velocity gradient between -35 and ~65 
km s" 1 kpc -1 . We adopted an intermediate value for the rotation- 
velocity gradient. 

In UGC 4483 the HI velocity dispersion <thi is only a fac- 
tor -2-3 smaller than the observed rotation velocity V rot . Thus, 
in order to trace the gravitational potential, the rotation curve 
has to be corrected for pressure support. We calculated the 
asymmetric-drift correction following Meurer et al. (1996). We 
assumed that the H I disk has constant scale-height and velo- 
city dispersion, and fitted the H I surface density profile (Fig. 6) 
with the Gaussian function Ehi (R) = £o X exp (-R 2 /2s 2 ^J, ob- 
taining So = 10.5 M Q pc -2 and s = 580 pc. The circular 
velocity V c i rc , corrected for asymmetric-drift, is thus given by 
^drc = Kot +(T Hi CK 2 /^ 2 )- Using the 3D models, we can constrain 
the mean velocity dispersion between ~6 and -10 km s" 1 . We 
assumed the intermediate value of 8 km s" 1 . The asymmetric- 
drift correction is significant only in the outer parts (see Fig. 4). 
Depending on the assumed value of cr m , the correction at the 
last measured points may be between ~2 and ~8 km s" 1 , giving 
a dynamical mass between about 1 and 2 xlO 8 M Q . 

In the above analysis, we used the centre and the P. A. derived 
from the stellar and HI morphologies. However, the differences 
between the observed and the best-model velocity fields (Fig. 5, 
top-right) show a systematic pattern that may point to a different 
location of the dynamical centre (cf. Warner et al. 1973, their 
Fig. 8). Thus, we built further models by changing the position 
of the centre and found that a shift of ~ 14" (-220 pc) towards the 
North- West reduces the pattern significantly. Moreover, a disk- 
model with an off-set centre and a P.A.^10° can reproduce the 
observed velocity field without radial motions (Fig. 5, bottom). 
However, these values for the centre and the P.A. are in marked 
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Fig. 6: HI surface density profiles (inclination corrected), de- 
rived separately for the northern (up-triangles) and southern 
(down-triangles) sides and for the entire galaxy (dots). 
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Fig. 7: R-band surface brightness profiles, derived separately for 
the northern (up-triangles) and southern (down-triangles) sides 
and for the entire galaxy (dots). The profiles are not corrected 
for internal extinction and inclination. The red-dotted line shows 
an exponential fit to the southern profile, giving yUo = 22.0 
mag arcsec" 2 (~ 40 Lo pc" 2 ) and# = 13.1" ± 0.1" (~ 203 pc). 



contrast with the observed stellar and HI morphologies. In par- 
ticular, the distribution of the old stars (which constitute an im- 
portant fraction of the dynamical mass in the inner parts, see 
Sect. 3.3) is remarkably symmetric and consistent with the cen- 
tre and P.A. previously adopted (see Fig. 6 of Dolphin et al. 2001 
and Fig. 9 of Izotov & Thuan 2002). Thus, in the following, we 
keep the dynamical centre coincident with the optical one. This 
has no significant effects on the derived rotation curve. 

3.3. Mass models 

In Sect. 3.2, we showed that UGC 4483 has a steeply-rising and 
flat rotation curve, indicating that there is a strong central con- 
centration of mass. To determine the relative contributions of 
luminous and dark matter to the gravitational potential, we built 
mass models following Begeman (1987). 
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The gravitational contribution of the stars was calculated us- 
ing a new approach: we estimated the stellar mass-to-light ratio 
M*/L using the HST results on the resolved stellar populations. 
Surface brightness profiles were derived from a sky- subtracted 
R-band image (Fig. 1, top-left), azimuthally averaging over a 
set of ellipses defined by (xo,yo), i, and RA (see table 1). The 
sky-background was determined by masking the sources in the 
field and fitting a 2D polynomial to the masked image. The 
surface brightness profiles were not corrected for internal ex- 
tinction, as UGC 4483 is extremely metal-poor (Skillman et al. 
1994) and the dust content is likely to be low, but were cor- 
rected for Galactic extinction assuming A R = 0.09 (Schlegel 
et al. 1998). Figure 7 shows the surface brightness profiles de- 
rived for the entire galaxy and for the southern (approaching) and 
northern (receding) halves separately. At R 20", the north- 
ern side of the galaxy is ~1 mag brighter than the southern 
one, owing to the recent starburst. Dolphin et al. (2001) and 
Izotov & Thuan (2002) resolved UGC 4483 into single stars 
and showed that the young stellar populations are mostly con- 
centrated to the North. The mass in young stars, however, con- 
stitutes only -10% of the total stellar mass: the SFH derived 
by McQuinn et al. (2010a) (modelling the color-magnitude dia- 
gram) implies that the stellar mass formed in the last 500 Myr 
is ~0.15xl0 7 M Q , whereas the stellar mass formed more than 
500 Myr ago is ~1.31xl0 7 M©. Thus, to calculate the gravi- 
tational contribution of the stars, we used the southern surface 
brightness profile, as it provides a better approximation of the 
stellar mass surface density. Assuming that 30% ± 10% of the 
gas is returned to the inter stellar medium by supernovae and 
stellar winds, the total stellar mass is (1.0 ± 0.3) x 10 7 M© 
(cf. McQuinn et al. 2010b), giving a stellar mass-to-light ratio 
M*/Lr = 1.1 + 0.3. The southern surface brightness profile can 
be fitted by an exponential-law with jjq = 22.0+0.1 mag arcsec -2 
(~ 40 L© pc" 2 ) and R = 13.1" + 0.1" (~ 203 pc). We assumed 
that the stars are located in a disk with vertical density distribu- 
tion given by p(z) = sech 2 (z/zo) (van der Kruit & Searle 1981) 
withzo 0.5R - 100 pc. 

The gravitational contribution of the gaseous disk was calcu- 
lated using the azimuthally-averaged H I surface density profile, 



which was derived from the total HI map. Figure 6 shows the 
H I density profiles derived by azimuthally-averaging over the 
entire galaxy and over the southern (approaching) and northern 
(receding) sides separately. There is a difference of the order of 
-2-3 M© pc" 2 between the two halves, that has a small effect 
on the resulting gravitational contribution. In agreement with the 
models in Sect. 3.2, we assumed an exponential vertical distribu- 
tion with zo = 100 pc. The mass of the HI disk was calculated 
from the total HI flux (see table 1) and multiplied by a factor 
of 1.32 to take into account the presence of helium. Molecular 
gas was not explicitly considered in the mass model because its 
amount is unknown. 

For the dark matter distribution, we assumed a pseudo- 
isothermal halo described by the equation 



where the central density po and the core radius r c are both free 
parameters of the mass model. 

Figure 8 (left) shows the rotation curve decomposition as- 
suming that M*/Lr=1.1, as is found by integrating the galaxy 
SFH (McQuinn et al. 2010a,b). The gravitational contribution 
of the stars V* is -50% of the circular velocity V c i rc at /? pea k = 
2.27?o - 450 pc. The parameters of the halo are: po = (101 + 20)x 
10" 3 M© pc" 3 and r c = 0.34 + 0.05 kpc. These values are com- 
parable with those found by S waters et al. (2011) for a sample 
of 18 dlrrs (assuming that M*/L R =1); the halo of UGC 4483, 
however, is one with the smallest core radius and highest cen- 
tral density. As discussed at the end of Sect 3.1, we applied the 
asymmetric-drift correction assuming a constant velocity disper- 
sion <thi = 8 km s" 1 . Different values of cr m would change the 
circular velocity at large radii and thus give slightly different va- 
lues for the dark matter halo. In particular, considering that <thi 
may be between 6 and 8 km s" 1 , the dynamical mass within the 
last measured point of the rotation curve may be between 1 and 
2 xlO 8 M Q and the baryon fraction between -43% and -22%, 
respectively. In any case, baryons (gas and stars) constitute a re- 
levant fraction of the total mass. 
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A stellar disk is defined to be maximum if F* = V*/Vcirc - 
0.85 ± 0.10 at 7? P eak (Sackett 1997; Bershady et al. 201 1). In our 
case F* = 0.5 ± 0. 1 ; the error takes into account the uncertainties 
on the stellar mass, the gas-recycling efficiency and the rotation 
velocity. This is in line with the results of the DiskMass survey 
(Bershady et al. 2011), that have measured the stellar velocity 
dispersion in a sample of spiral galaxies and found that stellar 
disks typically have F* ~ 0.5. We point out that the colour- 
magnitude diagrams of the resolved stellar populations provide 
the most direct method to quantify stellar masses, as the results 
depend only slightly on the assumed evolutionary tracks, me- 
tallicity, and sampling of the SFH (e.g. Annibali et al. 2003). 
The initial mass function (IMF), instead, may have a stronger 
effect: McQuinn et al. (2010a) assumed a single-slope Salpeter 
from 0.1 to 100 M© and a binary fraction of 35%. If the Chabrier 
and/or Kroupa IMFs are assumed, the stellar mass would sys- 
tematically decrease by, respectively, -25% and -30% and the 
gas-recycling efficiency may go up to -50%, further decreasing 
the stellar mass. 

We note that the molecular gas component has been ne- 
glected so far. UGC 4483 is undetected in the CO line with an 
upper limit of 0.195 K km s" 1 within the inner 55" (-800 pc) 
(Taylor et al. 1998). We extrapolated the relation between the 
CO-to-H2 conversion factor and metallicity from Boselli et al. 
(2002) down to the metallicity of UGC 4483 and obtained a 
corresponding upper limit on the H2 mass of Mh 2 < 10 7 M Q . 
Alternatively, we estimated the H2 mass from the empirical rela- 
tion between SFR and molecular mass, as determined by Leroy 
et al. (2008), although this relation may not apply to starburst 
galaxies. Assuming that SFR = 0.01 M© yr" 1 (McQuinn et al. 
2010a), we derived M H2 - 1.9 x 10 7 M©. Both estimates indi- 
cate that the molecular gas may be dynamically important and 
comparable in mass to the stellar component, implying a full 
maximum-disk situation, as we now discuss. 

Figure 8 (right) shows the maximum-disk solution (Sancisi 
& van Albada 1987). The stellar contribution can explain the in- 
ner parts of the rotation curve if M*/L R 3.3, i.e. if the stel- 
lar mass is -3 times higher than obtained from the study of 
the resolved stellar populations. This mass may be provided by 
molecules, if they are distributed in a way similar to the stars. 
We also find a good fit by fixing M*/Lr = 1.1 and scaling the 
HI contribution by a factor of -5, similarly to the results of 
Hoekstra et al. 2001 for other gas-rich galaxies. These results 
suggest that the distribution of the dynamical mass is closely 
coupled to that of the baryonic mass (see e.g. Sancisi 2004 and 
Swaters et al. 201 1). We also checked the position of UGC 4483 
on the baryonic Tully-Fisher relation and found that it follows 
the correlation within the observed scatter (cf. McGaugh 201 1). 

Finally, we considered the predictions of the Modified 
Newtonian Dynamics (MOND) (Milgrom 1983, see Famaey & 
McGaugh 2011 for a review). Using UGC 4483, we can test 
MOND without any free parameter, as the value of M*/L is 
provided by the colour-magnitude diagram of the resolved stel- 
lar populations and the distance is well determined from the tip 
of the red giant branch (Dolphin et al. 2001; Izotov & Thuan 
2002). We found that MOND systematically over-predicts the 
observed rotation curve by -5-6 km s" 1 . The discrepancy does 
not strongly depend on the assumed M*/L, as it is mainly driven 
by the atomic gas content, and it would further increase if 
molecules are also considered in the mass model. However, the 
discrepancy disappears if the inclination of the HI disk is as- 
sumed to be -43° instead of 58° (as is derived from the optical 
image). This possibility cannot be ruled out by the present ob- 
servations, because it is not possible to determine i neither from 
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Fig. 9: Comparison between UGC 4483 (blue line) and the dwarf 
irregular galaxy DDO 125 (UGC 7577) (red line). Top: HI sur- 
face density profiles. Middle: R-band surface brightness profiles, 
corrected for the inclination assuming that the disk is transpar- 
ent and using the ellipticity from Swaters & Balcells (2002). 
Bottom: HI rotation curves (not corrected for asymmetric drift). 

the total HI map (as the HI distribution is strongly lopsided) 
nor from the velocity field. 

4. Discussion 

In Sect. 3.2, we showed that UGC 4483 has a steeply -rising 
rotation curve that flattens in the outer parts at a velocity of 
-18-20 km s _1 . As far as we know, this is the lowest-mass 
galaxy with a differentially rotating H I disk. The steep rise of 
the rotation curve points to a strong central concentration of 
mass, which seems to be a characterizing property of BCDs. 
Indeed, similar results have been found also for other BCDs, 
e.g. I Zw 18 (Lelli et al. 2012), NGC 2537 (Matthews & Uson 

2008) , and NGC 1705 (Meurer et al. 1998). The steeply-rising 
and flat rotation curves of these BCDs are remarkable, as dlrrs 
typically have slowly-rising rotation curves (e.g. Swaters et al. 

2009) . The terms "steeply-rising" and "slowly-rising" are re- 
ferred to rotation- velocity gradients measured in physical units 
(km s" 1 kpc -1 ) that can be directly related to the dynamical mass 
surface densities in M pc -2 . If these gradients are expressed in 
terms of disk scale-lengths, rotation curves of dlrrs rise as steep 
as those of spiral galaxies (e.g. Swaters et al. 2009). 
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Fig. 10: Comparison between the properties of BCDs (squares) and of dlrrs (dots). The dlrrs are taken from Swaters et al. (2009), 
selecting the objects with high-quality rotation curves and V rot < 70 km s" 1 at the last measured point. The BCDs are UGC 4483 (this 
work), I Zw 18 (Lelli et al. 2012), NGC 6789 (Lelli et al. in preparation), and NGC 1705 (Meurer et al. 1998). Some galaxies are 
labelled according to their UGC number. Left: R-band central surface brightness (inclination corrected) versus the inner rotation- 
velocity gradient V rot CRo)/^o, where Ro is the optical scale-length. Right: peak HI surface density Shi, peak versus V rot (Ro)/Ro- 
Eh i, peak * s derived from the azimuthally-averaged H I surface density profiles. See Sect. 4 for details. 



In Lelli et al. (2012) (Fig. 10), we have illustrated the dyna- 
mical difference between the BCD-prototype I Zw 18 and a typi- 
cal dlrr of the same mass. Similarly, Fig. 9 compares UGC 4483 
with the dwarf irregular DDO 125 (UGC 7577, from Swaters 
et al. 2009). Assuming that DDO 125 is at a distance of 2.6 Mpc 
(Jacobs et al. 2009; Dalcanton et al. 2009), the HI observations 
of Swaters et al. (2009) have a linear resolution of -180 pc, si- 
milar to our linear resolution for UGC 4483 (-155 pc), making 
it possible to compare HI surface densities and velocity gra- 
dients. The two galaxies have approximately the same rotation 
velocity at the last measured point (-20 km s" 1 ) and dynamical 
mass (~l-2xl0 8 M©), but their structural properties are very dif- 
ferent: i) the H I surface densities of UGC 4483 are -3-4 times 
higher than those of DDO 125 in the inner regions (top); ii) the 
stellar component of UGC 4483 is much more compact than that 
of DDO 125 (middle); iii) the rotation curve of UGC 4483 has 
a steeper rise than the one of DDO 125 and flattens in the outer 
parts. These structural and dynamical differences between BCDs 
and typical dlrrs must be the key to understand the starburst phe- 
nomenon in BCDs. They also shed new light on the question of 
the progenitors and descendants of BCDs, as we now discuss. 

In Section 3.3, we showed that the central mass concentra- 
tion of UGC 4483 cannot be explained by the newly formed stars 
or by the concentration of H I. Old stars, although not sufficient 
to account for the inner rise of the rotation curve, constitute 
an important fraction of the mass in the inner parts. Similarly, 
the central mass concentration of I Zw 18 may be due to old 
stars and/or dark matter (unfortunately, the HST observations 
of I Zw 1 8 do not provide a direct estimate of the mass in old 
stars, see Contreras Ramos et al. 2011). For both galaxies, the 
gravitational contribution of molecules is very uncertain and it 
is unclear whether they are dynamically important or not. These 
results imply that either the progenitors of these BCDs are un- 



usually compact, gas-rich dwarfs, or there must be a mechanism 
that leads to a concentration of gas, old stars, and/or dark mat- 
ter, eventually causing the starburst. This mechanism may be ex- 
ternal (interactions/mergers) and/or internal (torques from mas- 
sive star-forming "clumps", see Elmegreen et al. 2012). It is also 
clear that, unless a redistribution of mass takes place, the descen- 
dants of BCDs must be compact dwarfs. Compact dlrrs do exist 
and have the following properties: i) high-surface-brightness ex- 
ponential profiles (see e.g. UGC 7690 and UGC 8550 in Swaters 
& Balcells 2002) or low-surface-brightness exponential pro- 
files with an inner light concentration (see e.g. UGC 6628 and 
UGC 12632 in Swaters & Balcells 2002), and ii) steeply-rising 
and flat rotation curves (see the same galaxies in Swaters et al. 
2009). The evolutionary links outlined here are in line with the 
results of e.g. Papaderos et al. (1996) and Gil de Paz & Madore 
(2005), which are based on surface photometry of large galaxy 
samples and indicate that the old stellar component of BCDs ge- 
nerally has higher central surface brightness and smaller scale- 
length than typical dlrrs and dSphs/dEs. We present here two 
plots that compare the structural and dynamical properties of 
BCDs and dlrrs. 

In order to compare the dynamics of BCDs and dlrrs, we use 
the inner rotation- velocity gradient V rot (R)/R as an estimate of 
central mass concentration. This is calculated at R = Ro, where 
Ro is the galaxy scale-length in the R-band; thus we take into 
account the different sizes of the systems. In Fig. 10, V rot CRo)/^o 
is plotted versus the central disk surface brightness fio (left) and 
the peak H I column density £Hi,peak (right) for a sample of both 
BCDs and dlrrs. The BCDs are UGC 4483, 1 Zw 18, NGC 1705, 
and NGC 6789. For I Zw 18 we used the R-band structural pa- 
rameters from Papaderos et al. (2002), the rotation curve from 
Lelli et al. (2012), and a distance of 18.2 Mpc (Aloisi et al. 
2007). For NGC 1705 we used the R-band structural parame- 
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ters from Gil de Paz & Madore (2005), the rotation curve from 
Meurer et al. (1998) named "model DD" (see their Fig. 10), 
and a distance of 5.1 Mpc (Tosi et al. 2001). For NGC 6789 
we used the R-band structural parameters and the rotation curve 
from Lelli et al. (in preparation), and a distance of 3.6 Mpc 
(Drozdovsky et al. 2001). The dlrrs were taken from the sam- 
ple of Swaters et al. (2009), by selecting the objects that have 
high-quality rotation curves (q < 2, see Swaters et al. 2009) 
and y rot < 70 km s" 1 at the last measured point. For Shi, peak, 
we used the peak value of the azimuthally-averaged H I density 
profile. The value of 2 H i, P eak depends on the linear resolution of 
the H I observations, but we checked that this does not strongly 
bias our analysis as the trend in Fig. 10 does not significantly 
change if one considers only galaxies with linear resolution < 
500 pc (corresponding to distances < 7 Mpc). We also point out 
that BCDs usually have clumpy H I distributions and locally the 
H I column densities can reach even higher values, up to 40-50 
M Q pc -2 , much higher than in dlrrs. 

The two diagrams in figure 10 show a clear trend: galaxies 
with a high rotation- velocity gradient have also high central sur- 
face brightness (left) and high peak HI surface density (right). 
The BCDs are in the upper parts of these distributions: they 
have V ro tCRo)/^o > 40 km s" 1 kpc -1 whereas dlrrs typically 
have V ro t(^o)/^o < 40 km s" 1 kpc -1 . However, some dlrrs 
(UGC 3966, UGC 8550, UGC 7232, and UGC 7690) have in- 
ner rotation- velocity gradients comparable to those of BCDs and 
thus, in this respect, are dynamically similar to BCDs. These ob- 
jects have peak HI column densities much smaller than those 
found in I Zw 18 and NGC 1705 (Fig. 10, right) and are can- 
didate progenitors/descendants of BCDs. It would be interesting 
to study the SFHs of these compact dlrrs and investigate if they 
experienced a starburst in the recent past. 

Finally, we have seen that the asymmetries in the velocity 
field of UGC 4483 can be described as a global radial mo- 
tion of 5 km s" 1 . Radial motions have been found in two other 
BCDs: NGC 2915 (Elson et al. 201 la,b) and I Zw 18 (Lelli et al. 
2012). Elson et al. (2011a) assumed that the spiral arms in the 
H I disk of NGC 2915 are trailing and concluded that the radial 
motions are an outflow. For I Zw 18 and UGC 4483, it is not 
possible to discriminate between inflow and outflow, as it is not 
known which side of the disk is the near one. For both galaxies, 
we calculated the timescales associated with the radial motions 
(- ^hi / V ra d) and found that they are of the same order of mag- 
nitude of the orbital times (^ ^Hi/^rot)- This suggests that any 
outflow or inflow must be very recent and possibly associated 
with the most recent burst of star-formation (~ 10-20 Myr, see 
McQuinn et al. 2010a for UGC 4483 and Aloisi et al. 1999 for 
I Zw 18). If the radial motions were an outflow, its kinetic en- 
ergy would correspond to only ~ 1 % of the energy released by 
supernovae. 

5. Conclusions 

We analysed archival H I observations of the blue compact dwarf 
galaxy UGC 4483 and built model datacubes to investigate its 
gas kinematics. Our main results can be summarized as follows: 

- UGC 4483 has a steeply-rising rotation curve that flattens 
in the outer parts at a velocity of -20 km s" 1 . This is, to 
our knowledge, the lowest-mass galaxy with a differentially 
rotating HI disk. Radial motions of ~5 km s" 1 may also be 
present. 

- The steep rise of the rotation curve indicates that there is a 
strong central concentration of mass. Mass models with a 



dark matter halo show that old stars contribute -50% of the 
observed rotation speed at 2.2 disk scale-lengths. Baryons 
(gas and stars) constitute an important fraction of the total 
dynamical mass. These conclusions are based on the stellar 
mass obtained from the color-magnitude diagram of the re- 
solved stellar populations. 
- The maximum-disk solution requires a stellar mass 3 times 
higher than observed, that could be provided by molecules. A 
good solution is also found by scaling the H I contribution by 
a factor of ~5. These results suggest that the distribution of 
the dynamical mass is closely coupled to that of the baryons. 

UGC 4483, together with other BCDs like I Zw 18 and 
NGC 1705, appears structurally different from typical dlrrs in 
terms of H I distribution, stellar distribution, and dynamics. In 
particular, a central concentration of mass (gas, stars, and dark 
matter) seems to be a characterizing property of BCDs. This im- 
plies that the starburst is closely related with the gravitational 
potential and the H I concentration. Our results also suggest that 
the progenitors/descendants of BCDs must be compact dwarf ga- 
laxies, unless a redistribution of mass (both luminous and dark) 
takes place before/after the starbursting phase. 
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